Structure and equation of state of interaction site models 
for disc-shaped lamellar colloids 
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We apply RISM (Reference Interaction Site Model) and PRISM (polymer-RISM) theories to 
calculate the site-site pair structure and the osmotic equation of state of suspensions of circular or 
hexagonal platelets (lamellar colloids) over a range of ratios of the particle diameter over thickness 
D/<j. Despite the neglect of edge effects, the simpler PRISM theory yields results in good agreement 
with the more elaborate RISM calculations, provided the correct form factor, characterizing the 
intramolecular structure of the platelets, is used. The RISM equation of state is sensitive to the 
number n of sites used to model the platelets, but saturates when the hard spheres, associated 
with the interaction sites, nearly touch; the limiting equation of state agrees reasonably well with 
available simulation data for all densities up to the isotropic-nematic transition. When properly 
scaled with the second virial coefficient, the equations of state of platelets with different aspect 
ratios D/a nearly collapse on a single master curve. 



I. INTRODUCTION 

Onsager's celebrated theory of suspensions of rod- 
like or disc-like cylindrical particles of diameter D and 
length (thickness) L predicts an isotropic to nematic 
phase transition, both in the rod limit L/D 3> 1 and in 
the plate limit L/D <C 1, on purely entropic grounds. 
While Onsager's prediction received abundant experi- 
mental confirmation for several lyotropic colloidal sys- 
tems of elongated rod-like particles, like the much studied 
TMV virus [31, the existence of such a transition in dis- 
persions of lamellar particles, first seen by Langmuir 1^, 
was only recently confirmed by experiments on polymer- 
grafted gibbsite platelets 3 . From a theoretical point of 
view, there exists no equivalent, in the case of platelets, of 
the Onsager limit, whereby the equation of state of hard 
rods is accurately described by the virial series truncated 
after second order in the limit L/D — > oo. The contri- 
bution of higher order virial coefficients to the equation 
of state of the cylinders remains large, even in the limit 
L / D — > (infinitely thin plates) , except at very low den- 
sities 0. For platelets of finite thickness [L/D ^ 0), 
virial coefficients of order higher than two {B^, B4, . . .) 
are not available. Hence there is a clear need to develop 
accurate theories for simple models of plate-like particles. 
Attempts have been made along various lines: 

— Scaled particle theory Q, which is very successful 
for hard sphere fluids, has been extended to prolate or 
oblate hard ellipsoids of revolution 0, with moderate 
success when gauged against Monte Carlo simulations Q . 
Much less work has been devoted to disc-shaped platelets 
and the results for infinitely thin platelets {L/D = 0) 3 
overestimate Monte Carlo data for the equation of state 
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quite significantly at densities approaching the isotropic 
to nematic transition [a. Scaled particle theory, or its 
extension by Boublik !lOl | , also overestimates the pressure 
of the closely related model of "cut hard spheres" , for 
which Monte Carlo data are available for several aspect 
ratios jHIH. 

— Onsager theory, based on the second virial coeffi- 
cient B2 alone, can be "rescaled" according to the pre- 
scription of Parsons fTsj and Lee Although this 
semi-empirical procedure gives reasonably good results 
for hard rod systems {L/D > 1), it is much less satisfac- 
tory for hard discs {L/D < 1) [l5l|. 

— Rosenfeld's very successful "fundamental measure 
theory" for hard sphere systems |0 has been extended 
to the Zwanzig model 0| for square platelets , but a 
direct comparison with disc-shaped particles is not pos- 
sible because within the Zwanzig model the platelets can 
only have three discrete orthogonal orientations. 

— Instead of modelling platelets as geometrically sim- 
ple convex bodies (like flat cylinders or cut spheres), an 
alternative is to consider platelets to be rigid, ordered 
arrays of interaction sites, distributed over the surface 
of a disc in some regular pattern [T9l |. Two examples 
considered in this paper are shown in Fig. ^ Site-site in- 
teractions are assumed to be spherically symmetric, i.e. 
to depend only on the distance between sites on differ- 
ent platelets. In Ref. the site-site pair potential was 
chosen to be of the screened Coulomb form, to mimic 
a dispersion of highly charged Laponite platelets j2Q|] in 
the presence of added salt, and the pair structure was 
determined by Molecular Dynamics simulations. 

In the present paper we return to this interaction site 
model, but replace the screened Coulomb interaction by a 
hard-sphere repulsion, in order to mimic as far as possible 
suspensions of uncharged, fiat, cylindrical platelets. The 
objective is to calculate the pair structure and osmotic 
equation of state from RISM (Reference Interaction Site 
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FIG. 1: Schematic representation of site patterns of model 
platelets investigated in this work. 



Model, m HI) and PRISM (polymer-RISM, [H 13) 
integral equations for the site-site pair distribution func- 
tions. A systematic comparison will be made between 
the results from the RISM and PRISM theories, in or- 
der to quantify the importance of the pre-averaging pro- 
cedure and of "edge" effects, which are involved within 
the simpler PRISM theory, where all sites are assumed 
to be equivalent ,2^]. Wherever possible, the theoreti- 
cal predictions will be compared to available simulation 
data. Parts of the present work were briefly reported 
elsewhere [25j . 

II. INTERACTION SITE MODELS 

The lamellar systems under consideration contain N 
disc-shaped platelets of diameter D — 2R inside a vol- 
ume V. Each platelet carries a regular pattern (or grid) 
of n interaction sites occupied by a hard sphere of diam- 
eter a. We consider the two grids pictured in Fig.^ The 
circular grid is made up of n,- concentric rings of diam- 
eter di = {i — 1){D — a)/{n,- — 1), with 1 < i < rir, 
comprising each 1 (for i — 1, i.e. at the origin) or 
6{i — 1) sites, adding up to a total of n — 7, 19, . . . , 
271 sites for = 2,3, ...,10. The grid has six-fold 
symmetry around the axis orthogonal to the disc, and 
along each of the in-plane symmetry axes, hard spheres 



on neighbouring sites overlap whenever the aspect ratio 
D/(j < 2(nr — 1). For small numbers of rings, it is clear 
that the resulting geometric object is very different from 
a convex body, since the interaction-site platelet will have 
holes, or a highly corrugated surface, depending on the 
relative values of D/a and n^. For any given D/a, the 
platelet constructed from concentric rings of spheres will 
tend to a cylinder of the same aspect ratio only in the 
limit rir oo. 

The hexagonal grid shown in Fig. ^ is identical to that 
or Ref. 19]. Concentric hexagons replace the rings, and 
the resulting total numbers of sites are identical to those 
obtained with the circular pattern. The polygonal shape 
is reminiscent of that of Gibbsite platelets used in recent 
exjj eriments at the Van't Hoff Laboratory in Utrecht 0, 

The interaction between two platelets is purely steric, 
and consists of the sum of hard sphere interactions be- 
tween all the pairs of sites on the platelets, say A and B: 

VAB ^"^^vilr.-Vjl) (1) 
ieAjeB 

where and are the positions of sites i and j, with 
I < i,j < n on platelets A and B, while v{r) is the 
familiar hard sphere "potential" : 

f \ / oo if r < cr , . 

= I if r > a 

This completely defines the interaction site model of 
platelets which we have investigated. The model is ob- 
viously athermal, and its reduced thermodynamic prop- 
erties will only depend on the platelet density p = N/V. 
The much-studied Laponite platelets have typical dimen- 
sions Z? ~ 25 nm and L (identified here with the diameter 
a of the sphere) c± 1 nm, corresponding to an aspect ra- 
tio D/(7 = 25 [231 ■ The results which will be presented 
here arc for aspect ratios 5 < D /a < 50, as well as for 
infinitely thin platelets (ct = 0). 



III. RISM AND PRISM THEORIES 

The pair structure of a fluid of identical particles, each 
carrying n distinct interaction sites, is characterized by a 
set of n(n -I- 1)/2 site-site intermolecular pair correlation 
functions hij{r) = gij{r) — 1. These are related to a set 
of n{n + l)/2 intermolecular direct correlation functions 
Czjir) by the RISM Ornstcin-Zernike (OZ) relations [2ll |. 
which are most conveniently written in matrix form and 
in Fourier space as: 

h(fc) = w(fc)c(fc)w(fc) + pw(fc)c(fc)h(fc) , (3) 

where h = [/ly(fc)], c = [cij{k)], and w arc n x n sym- 
metric matrices. The elements of w = [wij (fc)] are the 
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Fourier transforms of the intramolecular correlation func- 
tions. Provided the molecules are rigid, we have explic- 
itly: 

= (4) 

where is the bond length between sites i and j on 
the same particle. The RISM-OZ system Q must be 
supplemented by a set of closure relations relating the hij 
and Cij . Since only hard-core interactions are present, the 
simplest choice is the Percus-Yevick (PY) closure, which 
combines the exact core condition 

(^) = ~1 ; r <a (5) 

with the approximation 

Cy(r) = 0, r>a. (6) 

The number of independent hij and may be dras- 
tically reduced by obvious symmetry considerations. In 
the circular platelet model of Fig. ^ all sites on a given 
ring are equivalent, so that the total number of indepen- 
dent correlation functions is reduced from n{n + l)/2 to 
nr(nr -I- l)/2. Following Raineri and Stell zTj, we intro- 
duce the symmetry- reduced matrix W(/c) of intramolecur 
correlation functions with elements: 

where the indices / and J run over the number of in- 
equivalent classes of interaction sites and n/ or nj are 
the numbers of equivalent sites of classes / or J in the 
molecule. We also define the following direct correlation 
function matrix C with elements: 

Cjj{k) = nj ■ Cij{k) ■ nj (8) 

where Cij (fc) is any of the identical direct correlation func- 
tion of the original n x n matrix, with i G / and j G J. 
The reduced RISM-OZ equation now reads in riy x 
matrix form: 

h(fc) = W(fc)C(fc)W(fc) -I- pW(fc)C(fc)h(fc) , (9) 

which is once more supplemented by the PY closure re- 
lations (jSJ and ||HJ). In the case of the hexagonal grid 
in Fig. ^ the reduction procedure due to symmetry is 
slightly more elaborate. Consequently, we made the sim- 
ple (but slightly inaccurate) assumption that all sites on 
a given hexagonal shell (also referred to as "ring") are 
equivalent, which is strictly true only for the — 2 shell. 

Despite its diagrammatic inconsistency , RISM has 
proved to be a successful theory of the pair structure of 
many simple molecular fluids |22j . In the case of macro- 
molecular and colloidal systems, with very large num- 
bers of interaction sites, the number of coupled RISM 



equations become intractable, and a considerable simpli- 
fication follows from the assumption that all interaction 
sites are equivalent. This leads to the PRISM theory, 
first applied by Schweizer and Curro to long flexible poly- 
mers . PRISM theory neglects end effects in that case. 
In the case of rigid lamellar particles of finite diameter, 
the PRISM assumption is clearly less justified, partic- 
ularly so as the number n,- of concentric rings of sites 
increases. The set of correlation functions in Eq. ijJJ is 
replaced by a single, platelet-averaged function h, defined 
by: 

and similarily for c(fc). The platelet-averaged intramolec- 
ular correlation function is defined by: 

w{k) = ^^w,,{k), (11) 

and the resulting single PRISM-OZ relation reads 

h{k) = w{k)c{k)[w{k) + pnh{k)] , (12) 

where pn is the number density of interaction sites. Note 
that w{k) is proportional to the intramolecular form fac- 
tor F(k) of a single platelet: 

= l + ^T.'^^^ik)- (13) 

which characterizes the geometry of the distribution 
of sites, and hence the geometric shape of the model 
platelets. While the form factor accounts for the inter- 
ference of radiation scattered from different parts of the 
same particle in an X-ray or neutron diffraction experi- 
ment, the local order in the fluid or suspension is char- 
acterized by the set of pair distribution functions gij (r) , 
or, approximately, by the pre-averaged single combina- 
tion H12() . An important particular case is the centre-to- 
centre pair distribution function gcc('')- The experimen- 
tally accessible total structure factor is: 

S{k) = F{k) + ^ V %(fc) = F{k) + 5i„ter(fc) , (14) 

where the site-site partial structure factors Sij (fc) are the 
Fourier transforms of the gij (r) . The structure factor 
provides a direct link with thermodynamics, via the com- 
pressibility equation 29]: 

lim S{k) = pkBTxT , (15) 

where xt is the isothermal compressibility. The osmotic 
pressure then follows from: 

f^P^^= [ '^P [^(^' = 0)]-' . (16) 
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FIG. 2: Form factor of model platelets, Eq. II13L with 
D/g — 25 and different numbers n of sites, compared with 
the continuous distribution limit 1171 and the Lorentzian ap- 
proximant H18|l . Symbols are defined in the figure. 



Returning to PRISM theory, the OZ relation ifT^ 
must be supplemented by a closure relation for the pre- 
averaged functions h or c. In this work we use throughout 
the PY closure © and both within the RISM and 
PRISM representations. PRISM theory has so far been 
mostly applied to polymer solutions and melts, whence 
the name |2J|. For flexible polymers the form factor F{k) 
follows from a statistical average over single polymer con- 
formations, while in the case of rigid rods or platelets con- 
sidered in the present work, F(k), or equivalently w{k), 
are trivially known functions of the wave number k. The 
case of rigid rods was investigated within the PRISM 
formalism both as the limit of polymers with a persis- 
tence length exceeding the overall polymer length [30l |. 
and by using the correct form factor (see and refer- 
ences therein). 
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FIG. 3: Total structure factor for model platelets with D/a — 
25, and 91 (top) or 217 (bottom) sites. Circles: RISM; lines: 
PRISM results with the discrete [Eq. I^J, full Unes], or the 
continuous [Eq. 1171 . dashed lines in the bottom panel] form 
factors. Reduced densities are p* = pD'^ = 0.01, 0.40, 0.80, 
1.40, 2.80 (from top to bottom). 



IV. PAIR STRUCTURE 

We have solved the RISM and PRISM equations for 
various interaction site distributions, with number of 
sites n ranging from n = 37 {n^ — A) to n = 271 
(rir = 10), for aspect ratios D/cr ranging from 5 to 25 
and in the limit cr = 0. Test calculations with D/a = 
and tLt- — 18 have also been carried out. 

The advantage of using symmetry-reduced RISM and 
PRISM theories is that n may be increased within a 
reasonable computational effort, in order to mimic, as 
accurately as possible, a uniform distribution over the 
platelet area, and hence a smooth surface of a thin cylin- 
der. In fact, the limit of a continuous distribution may 
be achieved straightforwardly in the PRISM formalism. 

The importance of "holes" between sites on the grid 
may be seen in Fig. |21 where the form factors F{k), de- 



fined in Eq. are plotted for several values of the 

number of sites, and compared to the limit of a contin- 
uos distribution of sites, when the sum in Eq. H13|l goes 
over to an integral on the surface of the platelet, with the 
result |32|: 



F{k) 



{kRf 



1 



Ji{2kR) 
kR 



(17) 



where Ji denotes the cylindrical Bessel function of or- 
der one. The discrete nature of the site distribution is 
reflected in the irregular oscillations visible for ka > I 
in Fig. 121 In this regime, the correct fc^^ power law de- 
cay at large k is progressively recovered by increasing 
the number n of sites for a given diameter D. A model 
with rir = 9 (n = 217 sites), where the distance between 
nearest- neighbour sites is 1.5a, reproduces the continu- 
ous form factor H17|) with reasonable accuracy. Figure [21 
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FIG. 4: Intermolecular structure factor Sinter (fc), for model 
platelets with Dja = 25 or cx), calculated within PRISM, 
with the Lorentzian form factor [Eq. 1181 1. and the Bessel 



form factor [Eq. 117^ ]. Reduced densities are p* 
(upper curves) and 1.40 (lower curves). 



= 0.40 



also displays a simple Lorentzian approximation of the 
form factor, suitable for analytical calculation [25ll33l |: 



[kRY 



(18) 



Results for the overall structure factor (|14|l . and for 
the centre-to-centre structure factor obtained by solving 
the RISM and PRISM equations for the circular pat- 
tern of sites shown in Fig. ^ are presented in Figures 
13 For the hexagonal pattern, the full and approximate 
symmetry-reduced RISM calculations yield practically 
indistinguishable results (not shown), and the structure 
factors are very similar to those obtained with the circu- 
lar pattern, for the same number of sites n, aspect ratio 
D/a, and reduced density. Figure compares the RISM 
and PRISM results for the total structure factor at five 
different densities and for an aspect ratio D/a — 25. The 
same discrete form factor defined by Eqs. I^J and H13|) 
was used in the two integral equations for the models 
with n = 91 and 217 sites. In the n = 91 case, RISM 
and PRISM results are virtually indistinguishable, and 
this agreement is even better for smaller n. The agree- 
ment deteriorates somewhat as n increases, as illustrated 
by the case with n — 217, showing that the PRISM pre- 
averaging assumption gradually breaks down when the 
number of sites increases. As visible from Fig. 13 the 
PRISM structure factors are slightly modified by going 
to the continuum limit, characterized by the form fac- 
tor (|17|) . and move closer to the RISM data, pointing to 
some degree of cancellation of errors. 

A comparison of the intermolecular part of the struc- 
ture factors S'inter(fc) [see Eq. ()14f) ]. calculated within 
PRISM in the limit of a continuous distribution of sites 
{n — > 00) using the exact [Eq. (|17|l ]. and the approxi- 



mate [Eq. H18|l ]. form factors is made in Fig. 0] for size 
ratios a/D — and 0.04; the results are seen to be 
very sensitive to the choice of form factor, particularly 
at small wavenumbers fc, which will have implications for 
the equation of state calculated via Eqs. (|15|1 and (|16() . 

Note that the only justification for the use of the ap- 
proximate Lorentzian form factor is that it allows an 
analytic solution of the PRISM equation in the limit 
of infinitely thin {cr/D 0) hard as well as charged 
platelets '29|. Indeed, in that limit the PRISM equa- 
tion ({T^ . supplemented by the closure relation: 



c(r) = co(5(r) 



(19) 



can then be solved analytically [3^: the parameter cq is 
determined by enforcing the exact core condition h(y — 
0) = —1. The resulting structure factor reads: 



S{k) = w{k)+ph{k) 

= 2{2 + 2pA{R) + {kRfy^ 



(20) 



The result (PUll is, not surprisingly, identical to that ob- 
tained from the PRISM solution of the thread model, us- 
ing the same closure (I19II and the same approximant H18() 
to the form factor, with the radius R of the platelet re- 
placed by the radius of gyration i?g of the polymer [s^ . 
This coincidence illustrates the importance of using the 
correct form factor of the objects under consideration. 
The form factor F{k) for the platelets is known exactly 
[cfr Eq. ((T7|l ] , but only approximately in the case of poly- 
mers, where it is often taken to be of the Debye form, 
valid for ideal (gaussian) polymer coils js^. 

RISM results for the centre-to-centre structure factor 
Scc{k), calculated for the n = 217 sites model, are shown 
in Fig. 13 for two size ratios, D/a — 10 and 25, and 
four densities. As expected, the Scc{k) for the smaller 
size ratio show more structure because the corresponding 
effective packing fractions are larger. The correspond- 
ing centre-to-centre pair distribution function, g^dr), for 
D/a = 25, is shown in Fig.|Bl and compared to Eppenga 
and Frenkel's Monte Carlo (MC) data for infinitely thin 
platelets @ . The oscillations in the RISM results at short 
distances reflect the surface "roughness" due to the dis- 
crete site distribution, and are of course absent in the 
MC data. Bearing in mind the differences between the 
two models, the agreement between theory and simula- 
tion is reasonably good, and all curves exhibit the ex- 
pected cusp at r ~ £). At the highest density the MC 
data exhibit a steep rise followed by a striking inflection 
point around r ~ 0.4D, which is only faintly reproduced 
by RISM; for r < 0.4D, the MC gcc{r) lies well above 
the RISM prediction, pointing to some degree of local 
orientational ordering, which may be regarded as a pre- 
cursor to the isotropic-nematic transition taking place at 
pD^ ~ 4 0. Pair distribution functions gcdr) for the 
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FIG. 5: RISM centre-to-centre structure factors Scc{k), for 
platelet models with D/a = 25 (top) and 10 (bottom), for 
densities specified by the symbols. 



size ratio D/a ~ 10 are reported in the bottom panel 
of Fig. El for four densities. The quahtative behaviour 
is similar to that predicted by RISM for the larger size 
ratio; a similar cusp appears at r ~ D. 



V. THE EQUATION OF STATE 

From a knowledge of the site-site correlation func- 
tions, the only direct route to the osmotic equation of 
state is via Eqs. (fT^ and lfTE|l . Results from RISM 
and PRISM calculations for circular site patterns with 
37 < n < 271 sites, and an aspect ratio D/a = 25 are 
shown in Fig. [T] and compared to the MC data for in- 
finitely thin platelets 0. For any given reduced density 
p* = pD^, the pressure increases with the number n of 
sites. This is easily understood, because for lower values 
of n, the distance between neighbouring sites exceeds 2a, 
so that two platelets can, in fact, partially interpenetrate 
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FIG. 6: Top: RISM centre-to-centre radial distribution func- 
tion for platelets with D/a = 25 (lines, in order of increasing 
densities, from left to right). Symbols are the correspond- 
ing simulation data for infinitely thin hard discs Bot- 
tom: RISM radial distribution functions for platelets with 
D/a = 10. 



for favourable relative orientations. But for n = 217, 
the distance between neighbouring spheres (~ 1.5ct) is 
small enough to avoid interpenetration, and the value 
of the pressure saturates, as was checked by RISM test 
calculations with ten concentric rings, corresponding to 
n — 21\ sites. These saturated values of the pressure lie 
about 10% below the MC data in the intermediate den- 
sity range. Interestingly, the MC data show very little 
curvature of the equation of state for p* > 1.5, and fall 
below the predictions of RISM and of the fifth order virial 
expansion at higher densities. This "softening" of the 
equation of state may be indicative of local orientational 
ordering prior to the isotropic-nematic transition which 
occurs at p* ~ 4. Such local ordering, already apparent 
in the pair distribution function shown in Fig. El can- 
not be accounted for by RISM theory, which assumes full 
isotropy of the local correlations. The PRISM results, ob- 
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FIG. 7: Equation of state of hard platelets with D ja — 25 in 
the isotropic fluid range, as obtained through RISM (symbols) 
and PRISM (full lines) calculations with n = 37, 91, and 217 
sites and the same form factor of Eq. RISM test calcu- 

lations with n = 271, and PRISM results with the continuous 
Bessel form factor, Eq. 117L are also reported. Monte Carlo 
data and the 5th order viriai expansion for infinitely thin hard 
platelets are taken form Ref. Q- 



tained with the discrete form factor appropriate for 
a given number of sites, are in excellent agreement with 
the RISM pressures up to n = 91. Deviations between 
PRISM and RISM results become apparent for n = 127 
(not shown) and increasingly so for n — 217. However, if 
the continuum limit (n oo) is taken within PRISM, the 
form factor goes over to that given in Eq. (|17|l . and the 
resulting equation of state agrees then closely with the 
RISM results for n = 217, which are close to saturation, 
i.e. to the expected continuum limit within RISM. 

Similar RISM and PRISM results for the equation of 
state of discs with an aspect ratio D/a — 10 are com- 
pared in Fig. |H1 to MC data for "cut hard spheres" with 
the same aspect ratio 0. The trends are qualitatively 
the same as for the previous case with D/a = 25. The 
"softening" of the "exact" equation of state (as obtained 
by MC simulations) is again very striking at packing frac- 
tions T] > 0.2, where the MC data are bracketed by the 
predictions of the 5th order viriai expansion (which lie 
above) and the RISM predictions. Close to the expected 
isotropic-nematic transition at 77 ~ 0.3 11] the "satu- 
rated" RISM pressure (estimated from calculation for the 
n = 217 sites model) nearly coincides with the MC value. 
The continuum PRISM equation of state is close to the 
"saturated" RISM result as in the case with D/a = 25. 

The PRISM pressures for D/a = 5, 10, 15, 20, 25, 50 
and in the limit cr ^ are collected in Fig. O where the 
platelet density and the ratio P/k-oT are both reduced 
by the second viriai coefficient in the isotropic phase 0: 




FIG. 8: RISM and PRISM (with Bessel form factor) pre- 
dictions for the equation of state of model platelets with 
D/a = 10. Results are compared with Monte Carlo simu- 
lations of cut hard spheres of Ref. [Tl|. 
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FIG. 9: Equation of state of model platelets in reduced units 
of pressure and density, P* — j3PB2, p* = P-B2 (where B2 is 
the second viriai coefficient). Full line is a cubic fit of all data 
in the range 10 < D/a < (x>. 



B2 = ^Dy^{3 + n)D''a + ^Da\ (21) 

When P* — PPB2 is plotted versus p* — pB2, the data 
nearly collapse on a single "master curve" , which is rea- 
sonably well fitted by a cubic polynomial. Only the re- 
sults for D / a — h lie significantly above the master curve 
for p* > 1. 

Equations of state obtained with the circular and the 
hexagonal grids of Fig. ^ are reported in Fig. ^| for 
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FIG. 10: RISM equation of state for model platelets with 
different geometries, D/a = 25 and n — 217, 91, or 37 sites 
(from top to bottom). Full symbols; circular geometry; open 
symbols: hexagonal geometry (see Fig. 0. 77 = pvo is the 
dimensionless packing fraction and vq is the particle volume. 



FIG. 11: Equation of state of thin model platelets. RISM (full 
circles: Dja = 25, open circles: D/a = 50); PRISM with the 
Lorentzian form factor, Eq. (1181 (dashed lines, D/a = 25, 50 
and OQ, from top to bottom); PRISM with the exact Bessel 
form factor, Eq. 11711 (full lines, D/a = 25, 50 and 00, from 
top to bottom). 



VI. CONCLUSIONS 



D/a = 25, and n = 37, 91, and 217 sites. As visible, 
the pressures of both models agree closely, if the densi- 
ties p are properly reduced by the volume of a circular 
or hexagonal platelet in the continuum (n ^ 00) limit. 

We finally return to the prediction of the continuum 
PRISM theory, using the PY closure for finite thickness 
D/a = 25 and 50, and the closure (|19|l for infinitely 
thin (cr 0) platelets. PRISM results, obtained with 
the exact (|17f) and the approximate (|18|l form factors are 
reported in Fig. ^2 and compared to the RISM predic- 
tions for D/a = 2b and 50 (with nr = 9 and 18, respec- 
tively). Two main conclusions which may be drawn from 
the comparison are the following. First, the pressures are 
higher when a — 0.04Z3, compared to cr = 0.02£' (and 
cr = within PRISM), as expected for obvious excluded 
volume reasons, but the effect appears to be exagerated 
by PRISM theory. Finally, as already evident from the 
comparison of the PRISM structure factors shown in 
Fig. 21 the form factor has a considerable influence on 
the PRISM results with the Lorentzian approximant to 
the exact form factor always leading to an overestimation 
of the pressure. In view of the present discussion, the 
good agreement between the predictions of the contin- 
uum PRISM theory, obtained with the Lorentzian form 
factor H18jL and the "exact" MC data for infinitely thin 
platelets fsl noted in Ref. must be considered as 

partly accidental, due to a cancellation of errors arising 
from the combination of an inadequate form factor and 
the assumption of equivalence of all sites. 



We have presented results of what we believe to be the 
first systematic application of RISM and PRISM theo- 
ries to rigid, lamellar colloidal particles. While RISM 
has previously been almost exclusively applied to rigid 
molecules [l^l, and PRISM, as indicated by the first letter 
of the acronym, has been mostly applied to polymers |24| . 
our investigation presents a systematic comparison of the 
predictions of the two (related) theories as applied to 
colloidal platelets. In this paper we have restricted our 
attention to systems with excluded volume (hard core) 
interactions only between the sites distributed in a circu- 
lar or hexagonal pattern. For low numbers n of sites, the 
model platelets have many "holes" or exhibit highly cor- 
rugated surfaces, but for sufficiently large n, the spheres 
associated with neighbouring sites touch or overlap, and 
in the continuum limit n — s- cx), which can easily be dealt 
with within the PRISM theory, the topology of the multi- 
site plates goes over to smooth cylindrical or hexagonal 
bodies. 

The main findings and conclusions of the present work 
may be summarized as follows: 

(a) RISM theory predicts reasonable pair structure, 
as embodied in the centre-to-centre pair distribution 
function (cfr Fig. |SJ| or the total structure factor (cfr 
Fig. |2Jl , compared to available simulation data. The fully 
isotropic form factor used in all our calculations, and 
the ensuing complete decorrelation between plate orien- 
tations and the site-site correlation functions, do not al- 
low a proper description of the local orientational order- 
ing observed in the MC simulations. 

(b) The RISM equation of state converges rather slowly 
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to the hard platelet limit, as the number of interac- 
tion sites increases, but "saturation" is reached when 
the spheres associated with neighbouring sites are close 
enough to prevent the interpenetration of platelets. The 
limiting equation of state is strongly convex and shows 
no sign of the "softening" observed in MC simulations 
well before the isotropic-to-nematic transition is reached. 
The RISM pressures lie systematically below the predic- 
tions of a five-term virial series, and of the closely related 
Boublik equation of state. 

(c) The properly scaled RISM equation of state data 
over a wide range of aspect ratios D/cr collapse onto a 
single "master curve" . 

(d) The results of the much simpler PRISM theory, 
which assumes all interaction sites to be equivalent, agree 
unexpectedly well with their RISM counterparts, when 
the appropriate form factors are used. In other words, 
pre-averaging effects do not seem to play an important 
role. The continuum (n — » oo) limit of PRISM the- 
ory performs remarkably well, considering its simplicity 
and modest computational demands, provided the cor- 
responding exact form factor lfT7|l is used. The pres- 
sures obtained when the Lorentzian approximant (17) 
to the form factor is used in the PRISM calculations, 
are systematically too high. The substantial differences 
observed between the results of PRISM theory obtained 
with the exact and the Lorentzian form factor are reas- 
suring, since the latter has also been used for the thread 
model of polymers, and leads to the prediction that the 
equation of state of thin, non- intersecting polymers of ra- 



dius of gyration i?g would be the same as that of infinitely 
thin platelets of radius R — Rg. 

We are considering a number of extensions of the 
present work. In order to investigate the isotropic-to- 
nematic phase transition and other possible transitions 
to anisotropic phases, RISM and PRISM theories must 
be generalized to allow for anisotropic form factors, mir- 
roring the orientation of platelets relative to a direc- 
tor. The present work will also be extended to binary 
mixtures of platelets of different diameters (and hence 
numbers of sites) and aspect ratios. This will intro- 
duce a competition between depletion-induced demixing 
and orientational ordering, which is expected to lead to 
rich phase behaviour (15 , 1,8, 35] . A final extension of 
our work will be to charged platelets, modelled by us- 
ing interaction sites which carry charges [19]. This ex- 
tension is important for a description of aqueous disper- 
sions of natural, or synthetic clays, like the widely stud- 
ied Laponite, which exhibit poorly understood Coulomb- 
induced gelling behaviour [36j- Our preliminary investi- 
gations 18] encourage us to pursue a more systematic 
investigation of charged platelets within the framework 
of RISM and PRISM theory. 

ACKNOWLEDGEMENTS 

DC wishes to thank Prof. C. Caccamo for his con- 
tinuous encouragement and support during this work. 
Helpful discussions with Dr. R. Blaak are also gratefully 
acknowledged. 



[1] L. Onsager, Phys. Rev. 62, 558 (1942); Ann. N.Y. Acad. 

Sci. 51, 627 (1949). 
[2] For a review see S. Fraden in Observation, Prediction 

and Simulation of Phase Transitions in Complex Fluids, 

edited by M. Baus, L. F. Rull and J. P. Ryckaert (Kluwer, 

Dordrecht, 1994). 
[3] I. Langmuir, J. Chem. Phys. 6, 873 (1938). 
[4] F. M. van der Koij and H. N. W. Lekkerkerker, J. Phys. 

Chem. B 102, 7829 (1998). 
[5] R. Eppenga and D. Frenkel, Molec. Phys. 52, 1303 

(1984) . 

[6] H. Reiss, H. L. Frisch, and J. L. Lebowitz, J. Chem. Phys. 
31, 369 (1959). 

[7] For a review see M. A. Cotter in Molecular Physics of Liq- 
uid Crystals, edited by G. R. Luckhurst and G. W. Gray 
(Academic Press, London, 1979). 

[8] B. M. Mulder and D. Frenkel, Molec. Phys. 55, 119 

(1985) . 

[9] K. L. Savithramma and N. V. Madhusudana, Molec. 

Crys. Liq. Cryst. 74, 243 (1981). 
[10] T. Boublik, J. Chem. Phys. 63, 4084 (1975). 
[11] J. A. C. Veerman and D. Frenkel, Phys. Rev. A 45, 5632 

(1992). 

[12] S. D. Zhang, P. A. Reynolds, and J. S. Van Duijneveldt, 

Molec. Phys. 100, 3041 (2002). 
[13] J. D. Parsons, Phys. Rev. A 19, 1225 (1979). 
[14] S. D. Lee, J. Chem. Phys. 87, 4972 (1987); ibidem 89, 



7036 (1989). 

[15] H. H. Wensink, PhD dissertation, University of Utrecht 
(2004). 

[16] Y. Rosenfeld, Phys. Rev. Lett. 63, 980 (1989); Y. Rosen- 
feld, D. Levesque, and J. J. Weis, J. Chem. Phys. 92, 
6818 (1990). 

[17] R. Zwanzig, J. Chem. Phys. 39, 1714 (1963). 

[18] L. Harnau, D. G. Rowan, and J.-P. Hansen, J. Chem. 

Phys. 117, 11359 (2002). 
[19] S. Kutter, J.-P. Hansen, M. Sprik, and E. Boek, J. Chem. 

Phys. 112, 311 (2000). 
[20] Laporte Inorganics Laponite Technical Bulletin 

L104/90/A. 

[21] D. Chandler and H. C. Andersen, J. Chem. Phys. 57, 
1930 (1972); D. Chandler in Studies m Statistical Me- 
chanics, vol. 8, 275, edited by J. L. Lebowitz and 
E. W. MontroU (North Holland, Amsterdam, 1982). 

[22] For a review see P. A. Monson and G. P. Morriss, Adv. 
Chem. Phys. 77, 451 (1990). 

[23] K. S. Schweizer and J. G. Curro, Phys. Rev. Lett. 58, 246 
(1987); J. G. Curro and K. S. Schweizer, Macromolecules 
20, 1928 (1987). 

[24] For a review see K. S. Schweizer and J. G. Curro, Adv. 
Chem. Phys. 98, 1 (1997). 

[25] L. Harnau, J.-P. Hansen, and D. Costa, Europhys. Lett. 
53, 729 (2001). 

[26] F. M. van der Koij, E. Kassapidou, and 



10 



H. N. W. Lekkerkerker, Nature 406, 868 (2000); 
F. M. van dor Koij, and H. N. W. Lckkcrkorkcr, 
Phys. Rev. Lett. 84, 781 (2000); F. M. van dcr Koij, 
A. P. Philipse, and J. K. G. Dhont, Langmuir 16, 5317 

(2000) . 

[27] F. O. Raineri and G. Stell, J. Phys. Chem. B 105, 11880 

(2001) . 

[28] D. Chandler, R. Silbey, and B. Ladanyi, Molec. Phys. 46, 

1335 (1982). 

[29] J.-P. Hansen and I. R. McDonald, Theory of Simple Liq- 
uids, 2nd ed. (Academic Press, London, 1986). 

[30] G. T. Pickett and K. S. Schweizer J. Chem. Phys. 112, 
4881 (2000). 

[31] L. Harnau and P. Reineker, J. Chem. Phys. 112, 437 
(2000). 



[32] A. Guinior and G. Fournct, Small Angle Scattering of 

X-rays (Wiley, New York, 1955). 
[33] G. T. Pickett and K. S. Schweizer J. Chem. Phys. 110, 

6597 (1999). 

[34] See e.g. J. L. Barrat and J.-P. Hansen, Basic Concepts 
for Simple and Complex Liquids, p. 96 (Cambrisge Uni- 
versity Press, 2003). 

[35] H. H. Wensink, G. J. Vroogo, and H. N. W. Lekkerkerker, 
J. Phys. Chem. B 105, 10610 (2001); H. H. Wensink and 

G. J. Vroege, Phys. Rev. E 65, 031716 (2002). 

[36] A. Mourchid, A. Delville, J. Lambrad, E. Lecoher, 
and P. Levitz, Langmuir 11, 1942 (1995); H. Kellay, 

H. Tanaka, G. Wegdam, and J. Meunier, Langmuir 15, 
7534 (1999). 



